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Abstract 

Analytical and geometrical properties of generalized power-law 
(GPL) density profiles are investigated in detail. In particular, a one- 
to-one correspondence is found between mathematical parameters (a 
scaling radius, rg, a scaling density, po, and three exponents, a, /3, 
7), and geometrical parameters (the coordinates of the intersection 
of the asymptotes, xc, Vc-, and three vertical intercepts, 6, hp, b^, re- 
lated to the curve and the asymptotes, respectively): (ro, poi a, /3, 7) ^ 
{xc, Vc-, b, hp, b^). Then GPL density profiles are compared with simu- 
lated dark haloes (SDH) density profiles, and nonlinear least-absolute 
values and least-squares fits involving the above mentioned five pa- 
rameters (RFSM5 method) are prescribed. More specifically, the 
sum of absolute values or squares of absolute logarithmic residuals, 
Ri = log psDnifi) — log pgpl (?■«)) is evaluated on 10^ points mak- 
ing a 5-dimension hypergrid, through a few iterations. The size is 
progressively reduced around a fiducial minimum, and superpositions 
on nodes of earlier hypergrids are avoided. An application is made 
to a sample of 17 SDHs on the scale of cluster of galaxies, within 
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a flat ACDM cosmological model (Rasia et al. 2004). In dealing 
with the mean SDH density profile, a virial radius, rvir, averaged over 
the whole sample, is assigned, which allows the calculation of the re- 
maining parameters. Using a RFSM5 method provides a better fit 
with respect to other methods. The geometrical parameters, aver- 
aged over the whole sample of best fitting GPL density profiles, yield 
(a,/?, 7) (0.6,3.1,1.0), to be compared with (a,/?, 7) = (1,3,1), 
i.e. the NFW density profile (Navarro ct al. 1995, 1996, 1997); 
(a,/3,7) = (1.5,3,1.5) (Moore et al. 1998, 1999); (a,/3,7) = (1,2.5,1) 
(Rasia et al. 2004); and, in addition, 7 1.5 (Hiotelis 2003), deduced 
from the application of a RFSM5 method, but using a different def- 
inition of scaled radius, or concentration; and 7 ~ 1.2-1.3 deduced 
from more recent high-resolution simulations (Diemand et al. 2004; 
Reed et al. 2005). No evident correlation is found between SDH 
dynamical state (relaxed or merging) and asymptotic inner slope of 
the fitting logarithmic density profile or (for SDH comparable virial 
masses) scaled radius. Mean values and standard deviations of some 
parameters are calculated, and in particular the decimal logarithm of 
the scaled radius, ^vir, reads < log^vir >= 0.74 and (Ts log ^vir — 0-15- 
0.17, consistent with previous results related to NFW density profiles. 
It provides additional support to the idea, that NFW density profiles 
may be considered as a convenient way to parametrize SDH density 
profiles, without implying that it necessarily produces the best pos- 
sible fit (Bullock ct al. 2001). A certain degree of degeneracy is 
found in fitting GPL to SDH density profiles. If it is intrinsic to the 
RFSM5 method or it could be reduced by the next generation of high- 
resolution simulations, still remains an open question. 
keywords - cosmology: theory - dark matter - galaxies: clusters - 
galaxies: haloes. 

1 Introduction 

Recent observations of anisotropies in the cosmic microwave background put 
severe constraints on the cosmological parameters, by comparison with the 
predictions of current theories of structure formation and the evolution of 
the Universe. The addition of information from large-scale structure surveys, 
Hubble parameter determinations, and Type la supernova results, yields evi- 
dence for a fiat {^m + ^A ~ 1), low-density ~ 0.3; ~ 0.7) universe, a 
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Zeldovich power-law index of primordial fluctuations {us ~ 1), a (non bary- 
onic) dark matter density {Qdh^ ~ 0.16) dominant over baryon matter den- 
sity (Qfc/i^ ^ 0.023), a Hubble parameter (normalized to 100 km s~^Mpc~^) 
between one half and unity {h ^ 0.69), and a cosmological age between ten 
and flfteen billion years (T ^ 13.7 Gyr), which make the main information 
(Sicvcrs ct al. 2003). The above mentioned vahics arc consistent with the 
results deduced by different investigations (e.g., Rubino-Martin et al. 2003; 
Spergel et al. 2003; Fukugita & Peebles 2004). 

The related (ACDM) cosmological model (e.g., Diemand et al. 2004; 
Rccd ct al. 2005) is consistent with a bottom-up picture (hierarchical clus- 
tering) of dark matter haloes, where smaller systems formed first from initial 
density perturbations and then merged with each other to become larger sys- 
tems, or were tidally disrupted and accreted from previously formed larger 
systems. Most naturally, the density profiles of haloes are expected to be a 
two-parameters family. This is why, assuming that halo formation may be 
approximated to an acceptable extent by spherical collapse, each protohalo 
perturbation is characterized by two independent parameters e.g., mass and 
radius (or overdensity) , at some fiducial cosmological time. 

A successful two-parameter functional form for the halo profile, where a 
scaled density depends on a scaled radius, has been proposed by Navarro 
et al. (1995, 1996, 1997, the last quoted hereafter as NFW). They also 
argued that the density profile is universal, in the sense that its shape does 
not appreciably change (over two decades in radius) for different halo masses 
(spanning about four orders of magnitude) , initial density fluctuation spectra, 
or cosmological parameters. Many studies on the NFW "universal density 
proflle", both numerical (with increasing resolution) and analytical, were 
done after their proposal (for a more detailed discussion see e.g., Hess et al. 
1999; Klypin et al. 2001; Bullock et al. 2001; Fukushige & Makino 2001, 
2003; Miichet & Hoeft 2003; Zhao et al. 2003). 

Generally speaking, the NFW density proflle may be considered as a 
special case of the 5-parameter family (Hernquist 1990; Zhao 1996, 1997): 



related to the choice (a, (3, 7) = (1, 3, 1), where po and Tq are a scaling density 
and a scaling radius, respectively ^. 

^To the knowledge of the authors, the family of density profiles expressed by Eq. (1) 
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The density profile, expressed by Eq. (1), reduces to a power-law both 
towards the centre, r — > 0, and towards infinite, r — >• +00, where the expo- 
nent equals —7 and — /3, respectively. It may be conceived as a generalized 
power-law and, in the following, it shall be quoted as GPL density profile. 
On the other hand, matter distribution within a simulated dark matter halo, 
in the following, shall be quoted as SDH density profile. 

For fixed exponents, one among the two remaining free parameters, the 
scaling density and the scaling radius, may be related to the mass and the 
radius of the virialized region. For further details see e.g., NFW; Bullock et 
al. (2001); Rasia et al. (2004). 

Some doubts on the "universality" of the NFW density profile were cast 
by latter investigations. There are main orders of reasons against the idea of 
a universal, NFW density profile, which can be briefly summarized as follows. 

(i) A steeper slope in the central regions (e.g., Fukushige & Makino 1997, 

2001, 2003; Moore et al. 1998, 1999; Ghigna et al. 2000; but see 
Navarro et al. 2004, for a different point of view). 

(ii) A non universal slope in the central regions, which depends on the power 

spectrum of the initial density perturbation (Syer & White 1998), or on 
the mass (Jing & Suto 2000; Ricotti 2003). Additional support to this 
idea is provided by recent, high-resolution simulations (e.g., Fukushige 
et al. 2004; Navarro et al. 2004). 

(iii) A certain degree of degeneracy with regard to the exponents, {a, (3,^), 

in fitting various GPL to SDH density profiles, in the whole range of 
resolved scales (e.g., Klypin et al. 2001). 

(iv) Different criterions in fitting GPL to SDH density profiles, such as mini- 

mizing the maximum fractional deviations of the fit, max | log{pGPL/ Ph) — 
\og{psDH/ Ph)\ (e.g., Klypin ct al. 2001); the sum of the squares of 

was first defined by Hernquist [1990, Eq. (43) therein], even if liis attention was restricted 
to tlie special case, (a,/?, 7) = (1,4,1) whicli, in turn, was earlier proposed by Kuzmin 
& Veltmann (1973). A family of density profiles including the special case studied by 
Hernquist, the so called 7 models, where (a,/?, 7) — (1,4,7), "^^s given independently by 
Dehnen (1993) and Tremaine et al. (1994). A detailed investigation of general (a, /3, 7) 
models was performed by Zhao (1996, 1997). 
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absolute^logarithmic residuals, = E[log(pGPL/Pft) - ^og{psDH/ Ph)f 
(e.g., Bullock et al. 2001); the sum of squares of relative residuals, 
TX{psdh/ Ph- PgplI Ph)/{PGPL/ Ph)f (e.g., Fukushige & Makino 2003, 
2004), where ph is a normalization value. For a more detailed discussion 
see e.g., Tasitsiomi et al. (2004). 

(v) A gentler slope in the central regions (7 < 1) and a steeper slope suffi- 

ciently far from the centre (/9 > 3), under reasonable boundary condi- 
tions such as a finite halo mass and force-free halo centre, and a van- 
ishing density at infinite distance (Miicket & Hoeft 2003). Additional 
support to this idea is provided by recent high- resolution simulations 
(Navarro et al. 2004). 

(vi) A gentler slope in the central regions (7 ~ 1) to be consistent with ro- 

tation curves, deduced from the observations, of low surface brightness 
galaxies (McGaugh & de Blok 1988; de Blok et al. 2001), the Galaxy 
(Binney & Evans 2001), and dwarf galaxies (van den Bosh & Swaters 
2001). For a more detailed discussion, see Miicket & Hoeft (2003). 
In addition, the validity of the Jeans equation implies 1 < 7 < 3 for 
dark matter haloes (Hansen 2004), unless the effects of the baryonic 
component are taken into consideration (El-Zant et al. 2004; Hansen 
2004). 

(vii) A discrepancy with the dark matter distribution required to ensure hy- 
drostatic equilibrium of gas, deduced from measured X-ray brightness 
profiles, in clusters of galaxies (Arieli & Rephaeh 2003). 

Though the resolution of numerical simulations is increasingly high, still 
there is no general consensus, or insufficient investigation, about some ques- 
tions concerning dark matter halo density profiles, namely: (1) definition and 
formulation of universal density profiles (e.g., Huss et al. 1999; Bullock et al. 
2001); (2) connection between GPL and SDH density profiles; (3) dependence 
of GPL density profiles on the three exponents, (a, /5, 7) and the two scaling 
parameters, (ro,po); (4) extent to which two or more GPL density profiles 
fit the results of numerical simulations (e.g., Klypin et al. 2001; Fukushige 

^The term "absolute" here has not to be intended as "absolute value" , but as opposite 
to "relative". More precisely, y, — y{xi) is an absolute residual, while [yi — y{xi)]/y{xi) is 
the corresponding relative residual. 
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& Makino 2001, 2003); (5) degree of degeneration of the three exponents in 
fitting GPL to SDH density profiles (e.g., Klypin et al. 2001). 

The hierarchical collapse of dark matter into virialized haloes is hkely 
to have played a key role in the formation of large-scale objects, such as 
galaxies and clusters of galaxies. The halo profile has a direct dynamical 
role in determining the observable parameters of the baryonic subsystems. 
Therefore further investigation on the above raised questions appears to be 
important. 

To this aim, in fitting GPL to SDH density profiles, both a nonlinear 
least-absolute values and a nonlinear least-squares method are used in the 
current paper. The related boundary condition is that both the mass and 
the radius of the virialized region are determined by the computer outputs 
and the choice of the cosmological parameters. The main features of GPL 
and SDH density profiles arc outlined in sections 2 and 3, respectively. A 
comparison between GPL and SDH density profiles is performed in section 4. 
Nonlinear least-absolute values and least-squares fits are outlined in section 
5. The subject of section 6 is an application to a sample of 17 SDHs and the 
related mean density profile, on the scale of clusters of galaxies, taken from 
Rasiaetal. (2004). The results are then discussed. Some concluding remarks 
are drawn in section 7. Further investigation on a few special arguments is 
made in the Appendix. 



2 GPL density profiles 

Plotting GPL density profiles on a logarithmic plane, (0 log^ log/), neces- 
sarily implies use of dimensionless coordinates, defined as: 

e = - ; m) = - ; (2) 

where the scaled radius, ^, can be related to density profiles where the isopy- 
cnic surfaces i.e. constant density, are similar and similarly placed ellipsoids. 
For further details, see Caimmi & Marmo (2003). Accordingly, Eq. (1) re- 
duces to: ^ ^ 

independent of the scaling parameters. 
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The special choice: 

pt^2Vo ; rt = ro ; (4) 

translates Eq. (3) into: 

2x /? - 7 

•^(^) = ^7(l+^«)X ' ^ " ' 

which has an immediate interpretation, as: 

/(1) = 1 ; p(rt)=pt . (6) 

the scaling density, p\ coincides with the density on an isopycnic surface, 
where the radius equals the scaling radius, r — . For further details, see 
Caimmi & Marmo (2003). 

Scaled GPL density profiles, expressed by Eq. (3), include oo^ GPL den- 
sity profiles, expressed by Eq. (f), for the whole, allowed set of scaling pa- 
rameters, (ro,Po)- A similar situation occurs for polytropes (e.g., Caimmi 
1980). 

As Eqs. (2) and (3) imply null density at infinite radius, the mass distri- 
bution has necessarily to be ended at an assigned isopycnic surface, which 
defines a truncation radius, R. The mass within the truncation isopycnic 
surface is (Caimmi & Marmo 2003): 

47r 

M = M{R) - —rlpoVmas = MoUrnas ] (7) 

where Mq is a scaling mass and the profile factor, Umas, has the explicit 
expression: 

l^mas = ^ r fm'd^ ; (8) 

Jo 

and the integration is carried up to: 

= = - ; (9) 

which may be conceived as a scaled, truncation radius. 

The logarithmic GPL density profile, deduced from Eq. (3), is: 

iog/ = -7ioge-xiog(i + n ; (10) 
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it can be seen that the first and the third logarithmic derivative, calculated 
at log.^ = i.e. r = vq, yields: 

which discloses the geometrical meaning of the scaling radius. 

Geometrical meaning of the scaling radius in GPL density profiles 

With regard to logarithmic GPL density profiles, the maximum slope varia- 
tion rate occurs at the scaling radius, log^ = log(r/ro) = 0, where the slope 
equals the mean slope of the related asymptotes, —7 and —j3, respectively. 

In the special case of NFW density profiles, 7 = 1, /9 = 3, and the slope at 
the scaling radius equals —2 (e.g.. Bullock et al. 2001; Klypin et al. 2001; 
Hiotehs 2003). For a more detailed discussion, see Caimmi & Marmo (2003). 

In the limit of negligible values of the scaled radius, ^, with respect to 
unity, Eq. (10) reduces to: 

log/ = -7loge ; e«l ; (13) 

which represents, in the logarithmic plane (0 log ^ log /), a straight fine with 
slope equal to —7 and intercept equal to 0. 

In the limit of preponderant values of the scaled radius, ^, with respect 
to unity, Eq. (10) reduces to: 

iog/ = -/3ioge ; e»i ; (14) 

which represents, in the logarithmic plane (0 log ^ log /), a straight line with 
slope equal to — /3 and intercept equal to 0. 

It can easily be seen that the straight lines, expressed by Eqs. (13) and 
(14), meet at the origin and, in addition, represent the asymptotes of the 
logarithmic GPL density profile, expressed by Eq. (10). The special cases 
related to NFW and MOA (Moore et al. 1999) logarithmic density profiles, 
are plotted in Fig. 1. The above results hold for a > 0. The case o; < 
makes the asymptotes change one into the other. The limiting case a = 
makes either a vanishing density {(3 7^ 7) or the asymptotes coincide i.e. the 
curve reduces to a straight line (/3 = 7). 
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Figure 1: Logarithmic NFW (top) and MOA (bottom) density profiles (full 
curves), with their asymptotes (dashed hues), in the plane (Olog^ log/). 
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3 SDH density profiles 



Dark matter haloes simulations need three basic ingredients, namely: (i) a 
cosmological model with fixed values of the parameters; (ii) an environment 
with defined initial conditions; and (iii) an assigned computer code. The 
density profile during the evolution, is calculated through the following steps: 
(1) determine the centre of mass of the halo; (2) count the particles (bound to 
the halo) within spherical shells, centered on the centre of mass, and equally 
spaced in logarithmic distance i.e. log(ri+i/ri) =const; (3) evaluate the mean 
density within each shell. For further details see e.g., NFW; Klypin et al. 
(2001); Bullock et al. (2001); Fukushige & Makino (2001, 2003). 

Simulated haloes are characterized by a "virial" parameter, either the 
virial mass, Mvir, or the virial radius, Tyir, defined such that the mean density 
inside the virial radius is Ayir times the mean matter universal density, ph — 
Pcrit^m, at that redshift: 

4:71 2 

where Pcrit — 3if^/ (SttG) is the critical density for closure of the Universe. 
The critical overdensity at virialization, Ayir, is motivated by the spherical 
collapse model: it is below two hundreds for an Einstcin-dc Sitter cosmology, 
and exceeds three hundreds for a fiat ACDM cosmology where Urn ~ 0.3, at 
z^Q (e.g.. Bullock et al. 2001). 

Plotting SDH density profiles on a logarithmic plane, (0 log 77 log^), 
necessarily implies use of dimensionless coordinates, defined as: 

r p 

— ; i^iv) = — ; Ph = Pcru^m ; (i6) 

ryir Ph 

without loss of generality. 

An upper limit to the domain of SDH density profiles follows from the 
definition of virial radius: regions placed outside are still falling in, and their 
macroscopic kinetic energy has still to be converted into peculiar energy (e.g.. 
Cole & Lacey 1996; NFW). 

A lower limit to the domain of SDH density profiles is put by the oc- 
currence of numerical artifacts (mainly two-body relaxation) in the central 
regions, within about O.Olryir (e.g.. Bullock et al. 2001; Fukushige & Makino 
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2001, 2003, 2004; Navarro et al. 2004) or even less (e.g., Diemand et al. 2004; 
Reed et al. 2005). 

Accordingly, SDH density profiles appear to be closely related to the 
virialized region in the range: 

-2 < logr/ < ; (17) 

which shall be assumed in the following. A typical SDH density profile on 
the scale of cluster of galaxies, taken from a sample of 17 simulated haloes 
(Rasia et al. 2004), is represented in Fig. 2. Also plotted therein are the 
best linear fits, determined by use of a least-squares fit to simulated data, 
within the 7- region, —2 < log 77 < —1, and the /3-region, — 1 < log 77 < 0, 
respectively. For further details, see Appendix A. 



4 Comparison between GPL and SDH den- 
sity profiles 

A comparison between GPL and SDH density profiles necessarily implies that 
the truncation radius, and the mass enclosed within the truncation isopycnic 
surface, do coincide with the virial radius and the virial mass, R = Tyir and 
M — Myir, respectively. Then the combination of Eqs. (7) and (15) yields: 

Po^A^ ; (18) 

Ph ^mas 

Ur-— ; (19) 

where the scaled virial radius, .^vir, is usually defined as concentration in the 
special case of NFW density profiles (NFW). With regard to a generic GPL 
density profile, there are several definitions of concentration (e.g., Klypin et 
al. 2001). Throughout this paper we shall define the concentration as the 
scaled virial radius i.e. the ratio of the virial radius to the radius where 
the logarithmic slope of the density profile equals the mean slope of the two 
asymptotes, and the slope variation rate is maximum, according to Eqs. (11) 
and (12), respectively. 
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Figure 2: Logarithmic SDH density profile (open circles) for a typical dark 
matter halo on the scale of clusters of galaxies. The virialized region is safely 
represented in the range —2 < log 77 < 0. Also plotted are the best linear 
fits (dashed), determined by use of a least-squares fit to simulated data, 
within the 7-region, —2 < log 77 < —1, and the /^-region, —1 < log 77 < 0, 
respectively. 
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The comparison of scaled GPL density profiles, expressed by Eqs. (2), 
with scaled SDH density profiles, expressed by Eqs. (16), yields: 

^^Uv ; / = -^ ; (20) 

Po 

where ,^vir is defined by Eq. (19). Accordingly, a generic, scaled GPL density 
profile, expressed by Eq. (3), takes the equivalent form: 

and the related, logarithmic GPL density profile, is deduced by use of Eq. (18). 
The result is: 
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log = log Avir - log Umas + 3 log^v 

-7logevir-7logr;-xlog[l + (evirr;)"] ; (22) 

which depends on three exponents, (q;,/3, 7), and two scaling parameters, 
{tq, Pq). On the other hand, the scaled mass, I'mas, is defined by Eqs. (8) and 
(9); the virial radius, r^^^, is known from the computer run; and the critical 
overdensity, Ayir, is determined by the cosmological model. 

According to Eqs. (2), (11), (12), and (20), the maximum variation in 
slope occurs at r = ro i.e. log^ = i.e. log 77 = — log^vir- Then Eqs. (11) 
and (12) translate into: 

\ o ' / logr7=-log{vir 

where the slope at log 1] = — log ^vir equals the mean slope of the related 
asymptotes, —7 and —(3, respectively. 

In the limit of negligible values of the scaled radius, ^vir?7, with respect to 
unity, Eq. (22) reduces to: 

log i! = log A^ir - log Umas + 3 log ^vir 

-7log^vir - 7log?7 ; Cvirr7<l ; (25) 
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which represents, in the logarithmic plane (0 logrj log ■(/'), a straight line with 
slope equal to —7 and intercept equal to log Ayir — logUmas + (3 — 7) log^vir- 
In the limit of preponderant values of the scaled radius, ^vir?7, with respect 
to unity, Eq. (22) reduces to: 

log 1p ^ log Avir - log Umas + 3 log {vir 

-/31ogevir-/51ogr; ; {virr? > 1 ; (26) 

which represents, in the logarithmic plane (0 log r] log ip) , a straight line with 
slope equal to —j3 and intercept equal to log Ayir — log Umas + (3 — /3) log^vir- 
It can easily be seen that the straight lines, expressed by Eqs. (25) and 
(26), meet at the point {\ogri,\ogip) = [log(ro/rvii.), log(po/Pfe)], where Eqs. (18) 
and (2) have been used. In addition, the above mentioned straight lines rep- 
resent the asymptotes of the logarithmic GPL density profile, expressed by 
Eq. (22). Special cases related to NEW and MOA (Moore et al. 1999) loga- 
rithmic density profiles, are plotted in Fig. 3. 

5 The RFSM5 method 

Given a set of SDH density profiles, one is left with the problem of fitting 
a GPL density profile to each simulation and to the average on the whole 
set. To this respect, nonlinear fits shall be used, minimizing the sum of 
both absolute values and squares of absolute logarithmic residuals, Ri — 
^ogi/jsDHiVi) ~ log '0(3Pi(r7j), used in the literature (e.g., Klypin et al. 2001; 
Bullock et al. 2001). 

Strictly speaking, the problem reduces to a search of extremum points 
of minimum, with regard to a function, Y = F{Xi, X2, X3, X4, X^), which 
represents a 5-dimension hypersurface in a 6-dimension hyperspace. In gen- 
eral, no point of minimum can exist within the domain, if bounded, or a 
finite number, or infinite. If two or more minima are present, then degener- 
acy occurs. Owing to an intrinsic difficulty related to the above mentioned 
analytical procedure, a numerical alternative shall be exploited. 

More specifically, a 5-dimension hypergrid made of 10^ points is placed 
around a fiducial minimum, and the sum of both absolute values and squares 
of absolute logarithmic residuals, l-Ril and J2 Rh respectively, are evaluated 
at each point, and finally two (in general) distinct minima are localized. Then 
a new iteration is performed, with respect to a new hypergrid, centered near 
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Figure 3: Logarithmic NFW (top) and MO A (bottom) density profiles (full 
curves), with their asymptotes (dashed lines), in the plane (Ologr/ log'^). 
Values of ro/ryir and po/p/i have been arbitrarily chosen. 
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the minima, where the size has been reduced and superpositions on nodes 
of carher hypergrids have been avoided. For the calculations made in the 
current paper, three iterations have been revealed to be sufficient. 

In dealing with the hypergrid, it would be better to use parameters with 
an immediate geometrical meaning, instead of their analytical counterparts, 
{tq, pQ,a, f3,'-f). To this aim, the logarithmic GPL density profile, expressed 
by Eq. (22), has to be studied in detail, which is made in Appendix A. The 
geometrical parameters to be used, {xc,yc,b,bi3,b^), are the coordinates of 
the intersection of the asymptotes, and the vertical intercepts of the curve 
and the asymptotes, respectively. 

In summary, the procedure outlined above acts along the following steps. 

(i) Select a set of SDH density profiles, related to an assigned computer 
code and a specified cosmological model. 

(ii) Determine, for each SDH density profile, the GPL density profiles which 
minimize the sum of absolute values and/or squares of absolute loga- 
rithmic residuals, using a 5-dimension hypergrid in the 5-dimension 
hyperspace {OxcVcb bpbj). 

(iii) Perform the desired number of iterations around the minima, using 
a hypergrid with same number of points (but none in common), and 
reduced in size, in respect of its earlier analogon. 

(iv) Define a mean SDH density profile, over the whole set of simulations, 
and repeat the procedure used for a single SDH density profile. 

The result consists in a one-to-one correspondence between SDH density pro- 
files, including the related mean, and GPL density profiles which minimize 
the sum of absolute values and/or squares of absolute logarithmic residuals, 
in connection with the hypergrid used. In the following, the above mentioned 
procedure shall be quoted as RFSM5 (Residual Functions Sum Minimization 
within a 5-dimension hyperspace) method. The functions used in the cur- 
rent paper are the absolute value and the square of the absolute logarithmic 
residual. 
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6 An application to SDHs on the scale of 
cluster of galaxies 

Using a RFSM5 method, GPL density profiles are fitted to a sample of 17 
SDH density profiles, on the scale of cluster of galaxies within a fiat ACDM 
cosmology (Rasia et al. 2004, hereafter quoted as RTM). The values of the 
cosmological parameters used therein are: fl\ = 0.7; Qm = 0.3; Qh = 0.03; 
h = 0.7; as = 0.9; where the symbols have their usual captions (e.g., Klypin 
et al. 2001; Bullock et al. 2001) and, in particular, the indices m and 
b denote all matter (dark + baryonic) and baryonic, respectively. For a 
detailed discussion on the computer code, initial conditions, the resolution 
issues, and the way of finding the halo centre, see RTM and further references 
therein. 

Simulations include both dark and baryonic matter, but only the former 
is relevant to the aim of the current paper. Accordingly, the baryonic matter 
shall not be considered, and all the parameters shall be intended as related 
to the dark matter halo. 

The definition of the virialized region within each halo, via Eq. (15), needs 
the knowledge of the critical overdensity at virialization, Avir. With regard 
to total (dark + baryonic) matter, it depends on the cosmological model 
(e.g., Bullock et al. 2001) and, in the case under discussion, {Ayir)m = 323 
at z = 0, where all the sample haloes may be considered as virialized to an 
acceptable extent (RTM). 

If only the dark matter is considered, then (Avir)^ = A^ir = C(^?;ir)m; 
where C is the fraction of dark matter in each density perturbation, and 
averaging over the whole sample yields ( — 0.907 (RTM). Accordingly, the 
value: 

A™ = 0.907 ■ 323 = 292.961 ; (27) 

can be used to an acceptable extent ^. 

As clearly pointed out in RTM, owing to the random criterion used for 
selection, their sample haloes are characterized by varying dynamical propcr- 

^The above value of the critical overdensity was deduced from an earlier, unpublished 
version of RTM. It is slightly different from Avir = 0.903 • 323.7625 = 292.3576, deduced 
from the current version, which appeared when the calculations were performed in this 
paper. As the relative difference amounts to about 0.2%, the calculations were not repeated 
using the latter value of the critical overdensity. 
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ties: at the present time, some are more relaxed, while others are dynamically 
perturbed. The surrounding environment can also be quite different: some 
selected clusters are more isolated, while others are interacting with the sur- 
rounding cosmic web. Accordingly, the sample is expected to be good enough 
to provide unbiased conclusions, and the related modelling may be thought 
of as representative of an average cluster, in an average environment and 
dynamical configuration. 

6.1 Individual SDH density profiles 

The main features of sample haloes at 2; = are listed in Tab. 1. The mass, 
Mvir, has been taken from RTM, while the mass, M^j., has been deduced 
from Eq. (15). The apparent discrepancy between Mvir and M^r is owing 
to two different sources. First, a systematic contribution takes origin from 
the uncertainty on Avir and, second, a random contribution arises from the 
uncertainty on rvir, in both cases with regard to Eq. (15). An additional 
random contribution is related to averaging the fraction of dark matter over 
the whole sample, with regard to Eq. (27). 

The relative difference, \1 — M'^.^^/M^i^\, is less than one percent in all cases 
except 9, where it is less than one and half percent. Then the virial mass can 
be evaluated, to a good extent, by use of Eq. (15), taking the virial radius 
from the results of simulations. It is worth mentioning that the RFSM5 
method is independent of the value of the virial mass, while a change in 
the value of the virial radius makes SDH density profiles systematically shift 
along the horizontal direction, see Fig. 2. 

6.2 Averaged SDH density profiles 

Given a set of logarithmic SDH density profiles, the mean SDH density pro- 
file is obtained by averaging over the whole set the values related to each 
logarithmic radial bin, in the range of interest, expressed by Eq. (17). The 
value of the critical overdensity at virialization, Ayir, is fixed by Eq. (27), 
then a single free parameter remains: the virial radius, rvir, which allows the 
calculation of the virial mass, M^ir. The related values are determined by av- 
eraging over the whole sample, using the data listed in Tab. 1, and inserting 
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run 


type 


Af 

I V 


^vir 


M ■ 

iWyir 






no 


it 


ZoZO / ^ 


lyoo 


/ DOOU 


1 DU^U 


2 


5'02.10 


M 


278569 


2305 


125400 


125010 


3 


502.11 


M 


85159 


1553 


38340 


38234 


4 


503.05 


R 


294373 


2347 


132500 


131970 


5 


504.01 


R 


179681 


1991 


80820 


80565 


6 


504.07 


R 


146386 


1860 


65850 


65686 


7 


505.02 


M 


318653 


2197 


108700 


108249 


8 


506.01 


M 


427583 


2470 


153800 


153825 


9 


506.03 


M 


166855 


1796 


60020 


59136 


10 


507.01 


R 


275259 


1691 


49600 


49359 


11 


507.03 


R 


158345 


1407 


28530 


28433 


12 


508.01 


R 


190453 


1884 


68600 


68262 


13 


508.04 


R 


101482 


1529 


36560 


36489 


14 


509.03 


R 


159330 


1913 


71690 


71463 


15 


509.14 


R 


107229 


1675 


48250 


47971 


16 


510.03 


R 


58734 


1524 


36060 


36132 


17 


510.07 


R 


71937 


1628 


44170 


44045 



Table 1: Main features of sample haloes at 2; = 0. Column captions: 1 - case; 
2 - computer run; 3 - type {R - safely relaxed; M - safely a major merger 
occurring); 4 - number of dark matter particles within the virial radius; 5 - 
virial radius {h~^ kpc); 6 - virial mass (/i~^lO^°M0); 7 - virial mass deduced 
from Eq. (15). Both virial radii and virial masses are normalized to the 
dimensionless Hubble parameter at the current time, h. 
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the mean value of the virial radius into Eq. (15). The resuh is: 

Fvir = 1866 h-^ kpc ; Mvir = 66330 h'^ lO^^M© ; (28) 

and the application of a RFSM5 method yields best fitting GPL density 
profiles, with radius and mass equal to r^ir and M^^j., respectively. 

A mean virial radius has been preferred in place of a mean concentration 
(RTM) for the following reasons. First, virial radii are independent of GPL 
density profiles, contrary to concentrations, or velocity profiles, which should 
be calculated for any choice of the fitting profile. Second, the range of virial 
radii, 1407 < rvir/(/i~^kpc) < 2470, corresponds to relative errors of about 
25% and 33%, respectively, with regard to a mean value, fyir = 1866 h~^kpc. 
On the other hand, the range of concentrations (calculated for NFW density 
profiles), 5 < ^vir < 10, corresponds to relative errors of about 32% and 37%, 
respectively, with regard to a mean vahic, = 7.2976. Then the average 
of the virial radius should be preferred to this respect. 

Having in our hands a SDH density profile averaged over the sample, 
listed in Tab. 1, and a value of virial radius and virial mass, expressed by 
Eq. (28), we are left with the search of a best fitting GPL density profile. To 
this aim, six alternatives are exploited. The first one consists in the mere 
application of the RFSM5 method to the mean SDH density profile. 

Among the remaining five, three allow to fix the exponents in the GPL 
density profile, expressed by Eq. (1), and then minimize the sum of absolute 
values and/or squares of absolute logarithmic residuals, with respect to a 
single free parameter, the scaling radius, ro (e.g., Zhao ct al. 2003). Accord- 
ingly, the minimization is performed using a 2-dimension grid. The related 
procedure shall be quoted as RFSM2 (Residual Functions Sum Minimization 
within a 2-dimension space) method. The function used in the current paper 
is the absolute value and/or the square of the absolute logarithmic residual. 
The following special GPL density profiles are selected: NFW, MOA, and a 
best fitting profile deduced from both density and velocity distributions in 
sample haloes (RTM), hereafter quoted as RTM density profile "^.Thc related 
values of the exponents are (a,/?, 7) =(1, 3, 1), (3/2, 3, 3/2), (1, 5/2, 1), 
respectively. 

Finally, the last two alternatives among the six mentioned above, consist 
in calculating the mean values of the geometrical parameters, {xc, yc, b, bp, bj), 

Ht is worth mentioning that a different normalization has been used here for the scaUng 
density, po = {^virf^^{po)RTM- 



20 



over the whole sample of best fitting, GPL density profiles, with regard to 
the minimization of the sum of absolute values and squares of absolute log- 
arithmic residuals, respectively. 

6.3 Results 

As outhned in section 5, a RFSM5 method has been applied to each sample 
halo, hsted in Tab. 1, and to the related, averaged SDH density profile, which 
has been defined above. In addition, a RFSM2 method has been applied 
to the mean SDH density profile, in the special case of NFW, MOA, and 
RTM density profiles. The values of the exponents, a, x, Pi 1, the scaled 
radius, ^vir, the scaling radius, tq, the scaling density, po, and the sum of 
residual functions at the fiducial minimum, where f{Ri) — \Ri\ 

and f{Ri) — Ri, are listed in Tabs. 2 and 3, respectively. The following 
conclusions are deduced. 

(i) In general, different GPL density profiles best fit to an assigned SDH 
density profile, with regard to the minimization of the sum of absolute 
values or squares of absolute logarithmic residuals, respectively. The 
best fitting GPL density profiles in the two above mentioned alterna- 
tives, are closer each to the other in dealing with the mean geometrical 
parameters, than with the mean SDH density profile. 

(ii) The values of the exponents, (a,/?, 7), appearing in Eq. (1), and de- 
duced from the geometrical parameters averaged over the whole halo 
sample, {xc,yc,b,bp,bj), are {a,P,^) ^ (0.6,3.1,1.0). 

(iii) With regard to the mean SDH density profile, the GPL density pro- 
files which best minimize the sum of squares of absolute logarithmic 
residuals, occur in the following order of accuracy: (1) application of 
a RFSM5 method (ADP); (2) density profile related to mean values 
of geometrical parameters (AGP); (3) RTM; (4) MOA; (5) NFW. The 
best minimization of the sum of absolute values of absolute logarithmic 
residuals, imphes the same order as above, but with AGP and RTM 
density profiles interchanged. 

The SDH density profiles related to the current sample, listed in Tab. 1, 
and the related mean SDH density profile, are plotted in Fig. 4 (dots) together 
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case 
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7 




^^0 


lOVo 


E\R^\ 


1 


0.74154 


3.1652 


3.4698 


1.1226 


3.3884 


823.39 


0.51314 


0.37127 


2 


0.55959 


4.5587 


3.3698 


0.81879 


3.9446 


834.78 


1.6992 


0.67590 


3 


0.44184 


3.7405 


2.8307 


1.1780 


6.6069 


335.79 


2.6317 


0.77869 


4 


0.45956 


3.6948 


2.7388 


1.0409 


7.2577 


461.97 


2.6930 


0.76086 


5 


0.65993 


2.7518 


3.0723 


1.2563 


4.7315 


601.14 


0.74171 


0.45363 


6 


0.52739 


3.2392 


2.6766 


0.96825 


13.884 


191.38 


7.3923 


0.62044 


7 


0.59587 


4.0247 


3.4760 


1.0778 


4.6559 


674.11 


1.8207 


0.57012 


8 


0.66861 


3.9138 


3.4748 


0.85806 


4.1687 


846.45 


1.3190 


0.22446 


9 


0.56370 


3.9028 


3.1506 


0.95062 


6.4565 


397.38 


3.2377 


0.27807 


10 


0.42473 


3.8973 


2.7797 


1.1139 


6.1660 


391.78 


2.3455 


0.47658 


11 


0.65178 


3.7207 


3.5506 


1.1255 


3.2359 


621.15 


0.67645 


0.44895 


12 


0.66423 


3.5417 


2.9664 


0.61389 


7.9433 


338.83 


3.2752 


0.75673 


1 Q 

io 


U. /UzoU 


6. iZZo 


Q 1 OQ1 

o. izoi 






OOZ.OU 


i.ooUo 


U.Ooioy 


14 


0.74575 


3.2092 


3.2130 


0.81972 


5.1286 


532.86 


1.2029 


0.35644 


15 


0.82346 


2.1323 


2.7865 


1.0307 


6.5313 


366.37 


0.91250 


0.59466 


16 


0.71709 


3.8681 


3.7082 


0.93443 


4.0738 


534.43 


1.2889 


0.43596 


17 


0.47209 


3.1773 


2.6513 


1.1513 


5.7544 


404.16 


1.1488 


0.66217 


ADP 


0.54955 


4.5235 


3.3962 


0.91034 


3.8019 


701.18 


1.5496 


0.12966 


AGP 


0.60209 


3.4021 


3.0756 


1.0272 


5.5083 


483.96 


1.5422 


0.36819 


NFW 


1 


2 


3 


1 


6.35 


419.81 


0.90131 


0.85996 


MOA 


1.5 


1 


3 


1.5 


3.08 


865.52 


0.093933 


0.47639 


RTM 


1 


1.5 


2.5 


1 


13.050 


204.28 


2.1949 


0.32672 



Table 2: Parameters of GPL density profiles which (i) minimize the sum of 
absolute values of logarithmic absolute residuals, using a RFSM5 method 
with regard to 17 sample haloes listed in Tab. 1 (top), and (ii) fit to the 
mean SDH density profile, to a different extent (bottom). Cases correspond 
to computer runs in the former alternative, and to GPL density profiles in 
the latter. The GPL density profile which minimizes the sum of absolute 
values of absolute logarithmic residuals, using a RFSM5 method with regard 
to the mean SDH density profile, is denoted as ADP. The GPL density pro- 
file defined by geometrical parameters, {xc- yc- b, bg, b^), averaged over their 
counterparts listed on the top (first 17 rows), is denoted as AGP. The scaling 
radius, Tq, and the scaling density, po, are expressed in kpc and lO^°M0/kpc^, 
respectively. 
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A 





1 
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1 


85979 


2.7291 


3.4893 


1.1429 


3 6308 


768.43 


45733 


012843 


2 


70267 


2.6771 


2.8238 


0.94262 


7.1450 


460.86 


1.4462 


050984 


3 


50594 


3.6824 


3.0631 


1.2000 


4.4668 


496.68 


1.1755 


059798 

KJ * KJ \J t_y 1 t_y 


4 


0.40157 


4.8960 


2.9564 


99023 


4.1115 


815.48 


2.1391 


044834 

KJ * KJ ^L^LKJ ^-J 


5 


0.59295 


3.2455 


3.1244 


1.2000 


4.7863 


594.26 


1.0721 


0.015935 


6 


50437 


2.8516 


2.5152 


1.0769 


16.501 


161.03 


7.3632 


033055 


7 


0.60451 


4.1915 


3.5715 


1.0377 


4.2267 


742.56 


1.7070 


0.029137 


8 


0.66861 


3.9138 


3.4748 


0.85806 


4.1687 


846.45 


1.3190 


0.0049367 


9 


0.58477 


3.7717 


3.2083 


1.0027 


5.6234 


456.26 


2.2399 


0.0061089 


10 


0.43019 


3.7611 


2.7372 


1.1192 


6.6374 


363.95 


2.5132 


0.020731 


11 


0.69931 


3.3055 


3.4625 


1.1509 


3.7154 


541.00 


0.67645 


0.017674 


12 


0.67451 


3.0786 


2.7795 


0.70312 


9.1201 


363.19 


2.8856 


0.050487 


13 


0.68562 


3.3840 


3.1595 


0.83933 


6.8234 


320.12 


2.5718 


0.036570 


14 


0.81937 


2.4907 


2.9870 


0.94621 


5.3088 


514.77 


0.77666 


0.012273 


15 


0.70472 


2.8497 


2.8694 


0.86117 


7.6913 


311.11 


1.9963 


0.028821 


16 


0.74414 


3.9381 


3.8674 


0.93684 


3.7154 


585.99 


1.1226 


0.015259 


17 


0.44450 


4.0781 


2.8351 


1.0224 


4.6774 


497.23 


1.5144 


0.057718 


ADP 


0.56832 


4.0722 


3.3143 


1.0000 


3.6308 


734.22 


1.0238 


0.0014983 


AGP 


0.58866 


3.5008 


3.0528 


0.99204 


5.5170 


483.20 


1.5726 


0.0065805 


NFW 


1 


2 


3 


1 


6.44 


413.94 


0.93151 


0.046088 


MOA 


1.5 


1 


3 


1.5 


3.05 


874.03 


0.091827 


0.015492 


RTM 


1 


1.5 


2.5 


1 


13.350 


199.69 


2.3074 


0.0072902 



Table 3: Parameters of GPL density profiles which (i) minimize the sum of 
squares of logarithmic absolute residuals, using a RFSM5 method with regard 
to 17 sample haloes hsted in Tab. 1 (top), and (ii) fit to the mean SDH density 
profile, to a different extent (bottom). Cases correspond to computer runs in 
the former alternative, and to GPL density profiles in the latter. The GPL 
density profile which minimizes the sum of squares of absolute logarithmic 
residuals, using a RFSM5 method with regard to the mean SDH density 
profile, is denoted as ADP. The GPL density profile defined by geometrical 
parameters, {xc,yc,b,bg,b^), averaged over their counterparts listed on the 
top (first 17 rows), is denoted as AGP. The scaling radius, Tq, and the scaling 
density, are expressed in kpc and IO^^Mq/Icpc'^, respectively. 
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with their best fitting GPL counterpart (full curves), determined by use of 
a RFSM5 method in the range defined by Eq. (17). In most cases, GPL 
density profiles related to the minimization of the sum of absolute values 
and squares of absolute logarithmic residuals, are virtually indistinguishable, 
and sometimes coincident. 

Different fits to the mean SDH density profile, listed on the lower parts of 
Tabs. 2 and 3, are represented in Fig. 5, where log[(r/rvir)^(p/p/i)] has been 
plotted instead of \og{p/ph), to make different trends more evident. Curves 
related to the minimization of the sum of absolute values and squares of 
absolute logarithmic residuals are virtually indistinguishable with regard to 
GPL density profiles with fixed exponents: NFW, MOA, and RTM. 

The values of some analytical and geometrical parameters, tjadp-i related 
to the best fitting GPL density profile to the mean SDH density profile, are 
listed in Tabs. 4 and 5 together with their counterparts, rj, averaged over the 
best fitting GPL density profiles to the whole halo sample, via minimization 
of the sum of absolute values and squares of absolute logarithmic residuals, 
respectively. Also listed therein are the related standard deviations, as ,7, 
the standard deviations from the mean, Us rj, and the standard deviations 
from the standard deviation from the mean, as -ji, which are expressed as 
(e.g., Oliva & Terrasi 1976, Chap.V, §5.6.3): 
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(29) 

(30) 

(31) 
(32) 



logCvir, logi/^as, yC: ^^d. 



where n — 17; r] = a, P, 7, .^v 
owing to Eq. (39), log^vir = -xc- 

The following conclusions are deduced from Tabs. 4 and 5. 



(iv) Values of parameters, tiadp: related to the best fitting GPL density 
profile to the mean SDH density profile, are different from their coun- 
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Figure 4: The SDH density profiles related to the current sample, listed in 
Tab. 1, and the mean SDH density profile, denoted as ADP (open circles), 
together with their best fitting GPL counterparts (full curves), determined 
by use of a RFSM5 method in the range —2 < log(r/rvir) < 0. Two curves 
on each panel correspond to the minimization of the sum of absolute values 
and squares of absolute logarithmic residuals, respectively. The above men- 
tioned curves are virtually indistinguishable in most cases, and sometimes 
coincident. 
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Figure 5: Comparison between different fits (full curves) to the mean SDH 
density profiles (filled circles), listed on the lower parts of Tabs. 2 and 3. 
The function log[(r/rvir)^(p/ph)] has been plotted instead of \og{p/ph), to 
make different trends more evident. The vertical scale is related to the lower 
case. The remaining cases are, in turn, vertically shifted of one unity with 
respect to their closest counterpart, starting from the lower one, to gain 
clarity. With regard to case AGP, the upper and the lower curve correspond 
to the minimization of the sum of absolute values and squares of absolute 
logarithmic residuals, respectively, and the contrary holds for case ADP. 
Concerning the remaining cases NFW, MOA, and RTM (labelled here as 
RAS), the two curves are virtually indistinguishable. 
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3. 
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0. 


,31970 
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,077538 
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,013298 



Table 4: Comparison between parameters, tjadp, related to the best fitting 
GPL density profile to the mean SDH density profile, and their counterparts, 
rj, averaged over the best fitting GPL density profiles to the whole halo sam- 
ple, with regard to the minimization of the sum of absolute values of absolute 
logarithmic residuals. Also listed are the related standard deviations, as 
the standard deviations from the mean, cx^ ^, and the standard deviations 
from the standard deviation from the mean, as p. It is worth remembering 
that log^vir = —xc: according to Eq. (39). 
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.052454 


0, 


.0089958 



Table 5: Comparison between parameters, tjadp, related to the best fitting 
GPL density profile to the mean SDH density profile, and their counterparts, 
rj, averaged over the best fitting GPL density profiles to the whole halo 
sample, with regard to the minimization of the sum of squares of absolute 
logarithmic residuals. Also listed are the related standard deviations, as 
the standard deviations from the mean, ^, and the standard deviations 
from the standard deviation from the mean, as p. It is worth remembering 
that log^vir = —xc: accordiug to Eq. (39). 
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terparts averaged over the best fitting GPL density profiles to the whole 
halo sample, as expected from the theory of errors. 

(v) The exponents of best fitting, GPL density profiles, are close to their 

NFW counterparts, conform to [Nint(a), Nint(/3), Nint(7)] = (1,3,1), 
where Nint denotes the nearest integer. The difference increases from 
about one hundredth for 7 to about one tenth for /5, and to about one 
half for a. 

(vi) The rms error of the logarithm of the scaled radius, ^vir, is cr^ log^^;, = 
0.15 — 0.17, to be compared with Ug log^vir = 0.18 deduced from richer 
samples where (a,/?, 7) = (1,3,1); J = (0.5-1.0) xlO"/i-^m0; 11 < 
n < 14; and n is an integer (Bullock et al. 2001). 

6.4 Discussion 

The application of a RFSM5 method suceeds in minimizing the sum of abso- 
lute values and squares of absolute logarithmic residuals, with respect to GPL 
density profiles where the exponents are kept fixed, such as NFW, MOA, and 
RTM, which allows the following definition. 

Universal density profile. Let a RFSM5 method be applied to an as- 
signed set of SDH density profiles. The best fitting GPL density profile 
to the mean SDH density profile, is defined as the related universal density 
profile. 

In this view, "universal" has to be intended as nothing but "best fitting". 
Strictly speaking, the above statement should apply to a simulated halo 
sample which is representative of the whole set of real dark matter haloes. 

The minimization of the sum of absolute values or squares of absolute 
logarithmic residuals, makes a firm criterion for deciding which, among two 
or more GPL density profiles, best fits an assigned SDH density profile. The 
results of the current paper confirm earlier results about sample haloes on the 
scale of cluster of galaxies, namely (a) MOA density profiles provide a better 
fit with respect to NFW (Fukushige & Makino 2001, 2003), and (b) RTM 
density profiles provide a better fit with respect to NFW (RTM), with the 
additional result (c) RTM density profiles provide a better fit with respect 
to MOA. 
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A RFSM5 method has recently been used for determining the scaled 
radius (Hiotelis 2003), but using a different definition with respect to ^vir = 
fvu/fQ-i Eq. (19). In fact, the usual definition of concentration is c = rvir/r_2, 
where r_2 is the radius related to a logarithmic slope, defined by Eq. (49), 
dy/dx = -2 (e.g., Klypin et al. 2001; Bullock et al. 2001; Hiotelis 2003). 
The former definition seems to be more general, as it allows the maximum 
change in slope at the scaled radius, as pointed out in section 2. In addition, 
it makes the definition of concentration meaningful also in early times, where 
the slope of a GPL density profile may be smaller (in absolute value) than 2 
(Hiotehs 2003). 

In dealing with dark matter haloes on the scale of cluster of galaxies, 
Hiotelis (2003) finds GPL density profiles where the exponent, 7, attains a 
value of about 1.5, in contrast with the results of the current paper, j ^ 1. 
Such a discrepancy is probably owing to the different definitions of concen- 
tration, mentioned above, which have been used. 

On the other hand, values of asymptotic inner slopes of fitting density pro- 
files determined in the current paper, are consistent with their counterparts 
deduced from recent high-resolution simulations using a three-parameter fit 
involving scaling radius, scaling density, and asymptotic inner slope (6 sam- 
ple objects, Diemand et al. 2004) or a two-parameter fit involving scaling 
radius and asymptotic inner slope (16 sample objects. Reed et al. 2005). 
The related parameters are listed in Tab. 6, which shows agreement between 
different approaches, within the fiducial range, 7 =F ^cTs 7- All sample haloes 
are on the scale of cluster of galaxies, with the exception of three objects be- 
longing to the Reed et al. (2005) sample, which represent SDHs embedding 
the Milky Way and two dwarf galaxies, respectively. To get a homogeneous 
sample, a reduced set has been considered, which includes SDHs embedding 
only clusters of galaxies. A marginal discrepancy between results from the 
current paper and the original Reed et al. (2005) sample is found to disappear 
in dealing with the reduced sample, as shown in Tab. 6. 

With regard to the asymptotic inner slope of the logarithmic density pro- 
file, the current results conform to the validity of the Jeans equation, which 
demands 1 < 7 < 3 for dark matter haloes (Hansen 2004), but shallower 
slopes may occur if the effects of the baryonic component are considered 
(e.g., El-Zant et al. 2004; Hansen 2004). In addition, the inequality, 7 < 2, 
related to simple analytical treatments of dark matter haloes (Williams et al. 
2004), is also fulfilled by the current results, see Tabs. 2-5. On the contrary. 
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c 


n 


7 


(Tg J 


(Tg J 


A-7 


A+7 


AT7 


7 


7+ 


A 


17 


0.99973 


0.16188 


0.039261 


0.38584 


0.25657 


0.32120 


0.88195 


1.1175 


S 


17 


1.0018 


0.13893 


0.033696 


0.29868 


0.19820 


0.24844 


0.90071 


1.1029 


D 


6 


1.1617 


0.13732 


0.056060 


0.24167 


0.25833 


0.25000 


0.99352 


1.3299 


R 


16 


1.2875 


0.23910 


0.059774 


0.28750 


0.41250 


0.35000 


1.1082 


1.4668 


R 


13 


1.2692 


0.24285 


0.067353 


0.26923 


0.43077 


0.35000 


1.0671 


1.4713 



Table 6: Comparison between statistical parameters related to the asymp- 
totic inner slope, 7, deduced from different samples using different fits. Col- 
umn captions: 1 - case: A - current paper, minimization of the sum of abso- 
lute values of absolute logarithmic residuals; S - current paper, minimization 
of the sum of squares of absolute logarithmic residuals; D - Diemand et al. 
(2004); R - Reed et al. (2005); 2 - n: total number of sample objects; 3-7: 
arithmetic mean; 4 - cr^ ^: standard deviation; 5 - cx^ ^: standard deviation 
from the mean; 6 - A~7: maximum negative deviation from the mean; 7 - 
A"'"7: maximum positive deviation from the mean; 8 - A^7: mean maximum 
deviation from the mean; 9 - 7" = 7 — 3(Ts 7: lower limit assigned to the 
mean; 10 - 7'^ = 7 + 30-^ -: upper limit assigned to the mean. Sample haloes 
represent clusters or groups, with the exception of three objects belonging to 
the original R sample (excluded in the reduced R sample), which represent 
SDHs embedding the Milky Way and two dwarf galaxies, respectively. 
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the inequality, 7 < 1, related to force-free halo centre and vanishing density 
at infinite distance (Miicket & Hoeft 2003), is only marginally consistent with 
the current results. 

No evident correlation is found between SDH dynamical state (relaxed or 
merging) and asymptotic inner slope of the logarithmic density profile, —7, or 
(for SDH comparable virial masses) scaled radius, ^^ii., contrary to previous 
results (Ascasibar et al. 2004) related to a sample of 19 high-resolution SDHs 
on the scale of both clusters of galaxies (13 objects) and galaxies (6 objects), 
with regard to NFW density profiles. An investigation on richer samples 
could provide more information to this respect. 

The GPL density profiles which best fit to the averaged SDH density pro- 
file, are characterized by exponents, (a,/5, 7), satisfying [Nint(a), Nint(/5), 
Nint(7)] — (1,3,1), the last related to NFW density profile. But the corre- 
sponding deviations are not negligible with the exception of the one from 
7. The comparison with values averaged on the whole halo sample, discloses 
that the exponents, and 7, seem to fiuctuate around their NFW counter- 
parts, but the same does not hold for a, which has a mean value of about 
0.6. Accordingly, NFW density profiles cannot be conceived as universal, in 
the sense mentioned above, with regard to the current halo sample. 

On the other hand, following e.g.. Bullock ct al. (2001), NFW density 
profiles (or alternative functional forms) may be considered as a convenient 
way to parametrize SDH density profiles, without implying that it necessarily 
provides the best possible fit. This is why the scaled radius, ^vir, and the 
scaled mass, Vmasi can be interpreted as general structure parameters, not 
necessarily restricted to a specific density profile. In particular, any spread 
in ^vir and Vmas can be attributed to a real scatter in a "physical" scaling 
radius, defined by e.g., Eq. (11), rather than to inaccuracies in the assumed, 
"universal" density profile. For further details see e.g.. Bullock et al. (2001). 

Additional support to the above considerations is provided by the value 
calculated for the standard deviation of the decimal logarithm of the scaled 
radius, o-siog^vi^ = 0.15-0.17, which is very close to o"siog^^i^ = 0.18 deduced 
from a statistical sample of about five thousands of simulated haloes, within 
mass bins equal to (0.5-1.0) xlO"/i~^M0, where 11 < n < 14 and n is an 
integer (Bullock et al. 2001). 

The standard deviations listed in Tabs. 2 and 3, related to characteristic 
parameters of best fitting GPL density profiles, with regard to sample haloes, 
have been calculated under the assumption that they obey a Gaussian dis- 
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tribution, using Eqs. (29)-(32). The existence of a Gaussian distribution is 
a necessary, but not sufficient condition, for the vahdity of the central hmit 
theorem. In this view, the parameters under discussion are related to the 
final properties of the corresponding sample halo, which are connected with 
the initial conditions, ai, a2, oini intended as random variables, by a 
transformation, rf = a\ ■ ■ ■■■ ■ OLn- For further details, see Caimmi & 
Marmo (2004). 

In addition, it is worth of note that the application of a least-absolute 

vahics or a least-squares method (in particular RFSM5), in fitting GPL to 
SDH density proffies (e.g., Dubinski & Carlberg 1991; Klypin et al. 2001; 
Fukushige & Makino 2003) implies a (fiducial) Gaussian distribution of the 
SDH density profile, Dsdh = ^og{psDH/ Ph), around the expected value de- 
duced from the related GPL density profile, i/gpl — ^'^z{PgplI ph)-, for any 
fixed logarithmic radial bin centered on a; = log[(ri+i -|- rj)/(2rvir)]- It is the 
particularization, to the case of interest, of a well known result of the theory 
of errors (e.g., Taylor 2000, Chap. 8, §8.2). 

The results of the current paper confirm a certain degree of degeneracy in 
fitting GPL to SDH density profiles, as pointed out by Klypin et al. (2001). 
For instance, four GPL density profiles fit to the sample halo 5*02. 10, where 
the sum of absolute values and squares of absolute logarithmic residuals, 
and I]-Rf, the exponents, a, P, and 7, the scaled radius, ^vir, the 
scaled mass, I'mas, the scaling radius, Vq, and the scaling density, Pq, lie 
within the following ranges: 



0.675 <Y.\^^\ < 0-732 ; 

0.0509 <Y,Rf < 0.0537 ; 

0.559 <a< 0.727 ; 

2.73 <p < 3.37 ; 

0.818 < 7 < 0.943 ; 

3.89 < ^vir < 8.04 ; 

0.43K z/„„, < 3.10 ; 

409 < ro/kpc < 847 ; 

1.28 < 10Vo/(10^°Mo/kpc3) < 1.70 ; 



(33a) 
(33b) 
(33c) 
(33d) 
(33e) 
(33f) 
(33g) 
(33h) 
(33i) 



which could be explained in a twofold manner. 
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On one hand, a degeneracy could be intrinsic to the 6-dimension hy- 
perspace where the RFSM5 method works. Accordingly, the 5-dimension 
hypersurface, w — F{xc,yc,b,bi3,b^), defined by the sum of absolute val- 
ues or squares of absolute logarithmic residuals, happens to be parallel, in 
some finite region of the domain, to the principal 5-dimension hyperplane, 
(0 Xc ijc b bg which implies infinite extremum points of minimum^. 

On the other hand, a degeneracy could be owing to the restricted domain 
of SDH density profiles, defined by Eq. (17), which is limited by the virial 
radius on the right and by the occurrence of numerical artifacts (mainly two- 
body relaxation) on the left. In this view, a more extended range could reduce 
the degeneracy. As the fit must necessarily be restricted to the virialized re- 
gion, one shall wait for higher-resolution simulations involving five-parameter 
fits to test this possibility. But in recent high-resolution simulations two or 
three-parameter fits only have been used (e.g., Navarro et al. 2004; Diemand 
et al. 2004; Reed et al. 2005). 

The current attempt is limited to GPL density profiles, defined by Eq. (1), 
but different alternatives have been exploited in the literature, such as the 
Burkert (1995) density profile: 



which resembles the NEW density profile for r ~ 0.02rvir- The corresponding 
scaling and scaled radii may be related as: (ro)^ = (^o)ArFTy/1.52; (^vir)s = 
1.52(^vir)ArFW^; respectively (e.g.. Bullock et al. 2001). 

Another possibility is a profile that curves smoothly over to a constant 
density at very small radii (Navarro et al. 2004): 



where the parameter. A, prescribes how fast the density profile turns away 
from a power-law near the centre. In logarithmic form, Eq.(35) represents 

^An example in a 3-dimension space, (0.xyz), is provided by the surface of a cylinder: 
in the special case where the height, or a basis, is parallel to the principal plane, (Oxy), 
infinite extremum points of minimum occur. If otherwise, there is a single extremum point 
of minimum. 




(34) 




(35) 
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a class (defined by the parameter, — 2/A) of Sersic (1968) density profiles 
(Merritt et al. 2005). The best fit reads (19 sample objects, Navarro et al. 
2004): A T 3(7^ X = 0.17216 :f 0.021897. 

On the other hand, the RFSM5 method may be extended to any kind of 
density profile, keeping in mind that different classes could exhibit different 
geometrical properties. 

The "universality" of density profiles involving scaled parameters, has to 
be intended as in polytropes (e.g., Caimmi 1980). 

A single distribution in the abstract space of the scaled variables, = f{^), 
corresponds to oo^ distributions in the physical space, p/po = /(r/ro), pro- 
vided the free parameters are the scaling radius, ro, and the scaling density, 
Po, and the scaled radius reads ^ — t/tq. 

In deahng with SDH density profiles, it would be more germane to the mat- 
ter speaking about best fitting, instead of universal, density profiles. The 
validity of the fit has to be restricted to a fiducial range where simulations 
are not affected by spurious effects such as two-body relaxation, according 
to e.g., Eq. (17). In particular, the asymptotic inner slope is necessary for 
the definition of GPL density profiles, but any conclusion outside the above 
mentioned fiducial range may be at risk, in absence of some kind of (direct 
or indirect) observational support. 

7 Conclusion 

Analytical and geometrical properties of generalized power-law (GPL) den- 
sity profiles have been investigated in detail. In particular, a one-to-one 
correspondence has been found between mathematical parameters (a scaling 
radius, ro, a scaling density, po, and three exponents, a, {3, 7), and geo- 
metrical parameters (the coordinates of the intersection of the asymptotes, 
Xci Uci ^-iid three vertical intercepts, 6, 6^, 6^, related to the curve and the 
asymptotes, respectively): (ro, Po, a, (3, 7) ^ {xc, yc, b, bp, h^). 

Then GPL density profiles have been compared with simulated dark 
haloes (SDH) density profiles, and nonlinear least-absolute values and least- 
squares fits involving the above mentioned five parameters (RFSM5 method) 
have been prescribed. More specifically, the sum of absolute values or squares 
of absolute logarithmic residuals, i?j = ^ogpspH^'^'i) — log PGPL(^i), has been 
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evaluated on 10^ points making a 5-dimension hypergrid, through a few iter- 
ations. The size has progressively been reduced around a fiducial minimum, 
and superpositions on nodes of earlier hypergrids have been avoided. 

An application has been made to a sample of 17 SDHs on the scale of clus- 
ter of galaxies, within a flat ACDM cosmological model (Rasia et al. 2004). 
In dealing with the mean SDH density profile, a virial radius, r^ii., averaged 
over the whole sample, has been assigned, which has allowed the calculation 
of the remaining parameters. Using a RFSM5 method has provided a better 
fit with respect to other methods. 

The geometrical parameters, averaged over the whole sample of best fit- 
ting GPL density profiles, have yielded (a,/?, 7) ~ (0.6,3.1, 1.0), to be com- 
pared with (a,/5,7) = (1,3,1), i.e. the NFW density profile (Navarro et 
al 1995, 1996, 1997); (a,/?, 7) = (1.5,3,1.5) (Moore et al. 1998, 1999); 
(a,/9,7) = (1,2.5,1) (Rasia et al. 2004); and, in addition, 7 ^ 1.5 (Hiotelis 
2003), deduced from the application of a RFSM5 method, but using a dif- 
ferent definition of scaled radius, or concentration; 7 ^1.2-1.3 deduced from 
more recent fits (Diemand et al. 2004; Reed et al. 2005). No evident cor- 
relation has been found between SDH dynamical state (relaxed or merging) 
and asymptotic inner slope of the logarithmic density profile or (for SDH 
comparable virial masses) scaled radius. 

Mean values and standard deviations of some parameters have been cal- 
culated, and in particular the decimal logarithm of the scaled radius, ^vir, 
has been found to yield < log^vir >= 0.74 and CTsiogg^;, = 0.15-0.17, consis- 
tent with previous results related to NFW density profiles. It has provided 
additional support to the idea, that NFW density profiles may be considered 
as a convenient way to parametrize SDH density profiles, without implying 
that it necessarily produces the best possible fit (Bullock et al. 2001). 

A certain degree of degeneracy has been found in fitting GPL to SDH 
density profiles. If it is intrinsic to the RFSM5 method or it could be reduced 
by the next generation of high-resolution simulations, still remains an open 
question. 

Future work demands a generalization of the above results on two re- 
spects. First, the method could be applied to data related to dynamical 
mass distributions inferred in cluster of galaxies (or galaxies), and the fit 
be compared with its counterpart deduced from simulated density profiles. 
Second, the method could be applied using Sersic (1968) density profiles in 
place of the family defined by Eq. (1). An advantage is that both baryonic 
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and dark (non baryonic) mass distributions are well represented by the Sersic 
law (Merritt et al. 2005), which depends on four parameters instead of five. 
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9 Appendix 

9.1 A. Analytical and geometrical properties of loga- 
rithmic GPL density profiles 

The values of the vertical intercept related to the curve, b, and to the asymp- 
totes, 6/3 and b^, are readily determined by putting logrj = i.e. = 1 in 
Eqs. (22), (25), and (26), respectively. The result is: 

b = log Avir - log Umas + 3 log ^vir " 7 log Cvir " X log(l + Cvir) ; (36) 
by = log Avir - log Umas + 3 log ^vir " 7 log {vir ] (37) 
bfS = log Avir - log Umas + 3 log {vir - P log {vir ; (38) 

for sake of brevity, let us denote the intersection of the asymptotes in the 
logarithmic plane, (0 log?] logip), as C{xc,yc) where, owing to Eqs. (18), 
(19), (25), (26), (37), and (38), the explicit expression of the coordinates 
reads: 



which yields the following. 

Theorem. For a selected (but arbitrary) SDH density proGle, cosmological 
model, and GPL density profile, the intersection of the asymptotes in the 
logarithmic plane, (0 log 77 log'^), occurs at a point, C{xc,yc), where the 
coordinates are the decimal logarithm of the ratio between scaling radius and 
virial radius, and scaling density and mean (matter) density of the universe, 
respectively. 



xc^ - logCvir = - log — = log — ; 

'"0 '"vir 



(39) 




(40) 
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The combination of Eqs. (18), (19), (39), and (40) yields: 

log Vmas = log Ayir " Sxc " VC ] (41) 

accordingly, the vertical intercepts of the curve and the asymptotes, expressed 
by Eqs. (36), (37), and (38), reduce to: 

b = yc + ixc - X log [1 + expio(-Q;a;c)] ; (42) 
b^ = yc + Ixc ; (43) 
= 2/c + ; (44) 

where, in general, exp^x = , and expx = e^, according to the standard 
notation. 

A change of coordinates, defined as: 

x = log77 ; y = logV' ; (45) 

makes the expressions of the curve and the asymptotes, defined by Eqs. (22), 
(25), and (26), translate into: 

y = yc-l{x-xc)-x^og{l + eyi.^^Q[a{x-xc)]} ; (46) 
y = yc - - xc) ; X <^xc ; (47) 
y = yc - f3{x - xc) ; x > ; (48) 

with regard to the plane (Oxy). 

The first derivative of the curve, defined by Eq. (46), takes the expression: 

^ = + ^ ; (49) 

dx 1 + exp^Q[Q;(a; — Xc)] ' 

and the particularization of the above result to the vertical intercept of the 
curve, reads: 

{¥\ =-/5 + — V ; (50) 

ydx/^^Q 1 + expio(-Q;a;c) 

according to Eqs. (42) and (50), the equation of the tangent to the curve at 
the vertical intercept is: 



y = 



expio(-aa;c) 



x + 6 ; (51) 



41 



the intersections of this hne with the asymptotes, expressed by Eqs. (47) and 
(48), let they be M^xmtVm) and R{xR,yR), respectively, may be calculated 
after some algebra. The result is: 

xm -| [1 + expio(-Q;xc)] ; (52) 

P - 7 

Vm ^bfj- (3 1^ _ [1 + ex.p^o{-axc)] ; (53) 
b -b 

xr = [1 + exp^o{-axc)] ; (54) 

P - 7 

yR^b^ + 7-^^ [1 + expiQ{-axc)] ; (55) 

P - 7 

the curve, the tangent to the curve at the vertical intercept, and the asymp- 
totes, are represented in Fig. 6. It is apparent that the curve must necessarily 
lie below the asymptotes and the segment, MR, with the exception of the tan- 
gential point, B(0,6). 

The combination of Eqs. (42), (43), and (44) yields: 

^7 - ^ = X log [1 + expio(-aa;c)] ; (56) 
b0-bj = {(3- 'y)xc = x^xc ; (57) 

and the comparison between the alternative expressions of the exponent, x, 
deduced from Eqs. (56), (57), produces: 



Aw = log [1 -I- expj^o ''^] (58a) 
A = > 1 ; (58b) 

w — —axc > ; (58c) 



which is equivalent to: 



u^-u- 1^0 ; (59a) 
u — exp^Q w ; (59b) 

where A > 1, w > 0, and — /3 < —7 < 0, in the case of interest (logarithmic 
GPL density profiles of the kind represented in Fig. 6). 
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Figure 6: A NFW logarithmic density profile, together with the tangent at 
the vertical intercept and the asymptotes, on the {Oxy) plane. The above 
mentioned straight lines define a triangle, CMR. The curve must necessarily 
lie below the asymptotes and the segment, MR, with the exception of the 
tangential point, B. 
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The function, 4>{u), on the left-hand side of Eq. (59a), has the following 
properties: 



(p{u) = u^ -u-l ; 0<M< +00 ; (60a) 
0(0) = -1 ; = ; (60b) 

lim (j){u) — +00 Muq) — ; (60c) 

u— ♦+00 

where Umin and Uq denote the abscissa of the extremum point (of minimum) 
and the zero of the function, respectively. Then Eq. (59a) may be solved by 
use of an iterative method. 

The combination of Eqs. (57), (58c), and (59b) yields: 

a — ; (61) 

xc 

X = T — ; ^2) 

\oguo 

on the other hand, the exponents, 7 and /3, may be deduced from Eqs. (43) 
and (44), respectively, as: 

(5=^-^^ ; (64) 

and the scaling parameters, ro and po, may be deduced from Eqs. (39) and 
(40), respectively, as: 

= r^ir expio xc ; (65) 
Po = Ph expio yc ; (66) 

the set of Eqs. (61)-(66) yields the following. 

Theorem. For a selected (but arbitrary) SPH density profile, cosmological 
model, and GPL density profile, in the logarithmic plane, (Oxy), there is 

a one-to-one correspondence between analytical parameters, {tq, p^.a, f3,'~f), 
and geometrical parameters, {xc,ycib,bp,h^), in the sense that either set 
univocally determines a GPL density profile. 
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The advantage of using geometrical instead of analytical parameters lies in 
a better understanding of how the curve change as one or more parameters 
do. 

With regard to SDH density profiles, according to Eqs. (17) and (45), let 
us divide the domain into two distinct regions, as: 

-2<x^<-l; -l<xp<0 ; (67) 

which shall be called the 7 region and the /3 region, respectively. 

Let (xj, yi) be coordinates of a generic point of a logarithmic SDH density 
profile, and [xi, y{xi)] their counterparts related to a fitting, logarithmic GPL 
density profile. Owing to Eq. (46), the corresponding, logarithmic absolute 
residual, is: 

Ri = yi- y{xi) =yi-yc + - xc) 

+ xlog{l + expio[a(a;i - xc)]} ; (68) 

the particularization of Eq. (68) to the 7 and /? region, defined by Eq. (67), 
allows the application of a least-squares fit to the related portions of SDH 
density profile. The best hnear fites are: 

y = bs-r - IsX ; (69) 

y = bsf3- PsX ; (70) 
and the coordinates of their intersection point, Cs{xs,ys), are: 

bs-r — bsfi bs^Ps — bsjsls \ 

= -7> ; Vs = — n ; 71 

Ps -Is Ps- Is 

the best linear fits to a selected SDH density profile, are plotted in Fig. 2. 

The intercepts, 6^7 and hsf}-, and the slopes, —75 and —/3s, appearing 
in Eqs. (69) and (70), arc calculated using the standard formulation of the 
method (e.g., Secco 2001, Chap. 4, §4.1): 



-A. = ; (72) 

X\Xx - XxXx 



, yxXxXx - XxXxXx 

bsX = ^ — ; (73) 

XxXx - XxXx 

XxXx - XxXx 



^2 ^2 XxXx /_ s 

XxXx -XxXx 
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where A = 7, a bar denotes arithmetic mean over the corresponding range 
of simulated values, and is the empirical variance of the histogram calcu- 
lated for the selected random variable. 

The comparison between best linear fits to SDH density profiles, defined 
by Eqs. (69) and (70), and related asymptotes of GPL density profiles, de- 
fined by Eqs. (47) and (48), implies the following, fiducial conclusions: (i) 
SDH best linear fits lie below related GPL asymptotes, in the range of in- 
terest, and (ii) SDH best linear fits are more inclined (in absolute value) 
with respect to related GPL asymptotes towards negative infinite, and less 
inclined (in absolute value) with respect to related GPL asymptotes towards 
positive infinite. Accordingly (iii) the intersection between SDH best linear 
fits lies below, and on the left, with respect to the intersection of related GPL 
asymptotes. 

The above mentioned conclusions (i) and (ii) read: 

b <bs(} ; bs/3 < bp ; b^^ < b^ ; (76) 
7, < 7 < ; (3</3s<-fs ; (77) 

on the other hand, GPL asymptotes intersect within the ( — h) quadrant, 
according to Eqs. (39) and (40), which yields: 

b<yc ] bp <yc ; < yc ; (78) 

finally, lower and upper values to the ranges: 

and the lower value to the range: 

bmin <b<bs/3 ; (80) 

may be deduced from SDH density profiles. The upper value of inequality 
(80) is owing to the first inequahty (76). The combination of inequahties 
(76) and (78) yields: 

bsp <b/3 <yc ; bs-y <b^ <yc ; (81) 
due to negative slopes of GPL asymptotes, —7 < and — /3 < 0. 
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With regard to a 5-dimension hyperspace, (0 xc yc b ^7)1 inequalities 
(79), (80), and (81), define a 5-dimension hyperrectangle of sides {xmax — 

^min)j {ymax~ Vminj^ (ps()~bmin)^ {Vmax ~ bgp), and {Umax ~ bgj), respectively. 

Then it is possible to make a 5-dimension hypergrid where the points are 
equally spaced; according to inequalities (81), the ranges yc < bp < ymax 
and yc < < ymax are to be excluded as outside the cases of interest. 
In dealing with the remaining points, the sum of both absolute values and 
squares of absolute logarithmic residuals, defined by Eq. (68), is performed, 
and the minimum on the hypergrid is determined in each alternative. Let it 
be (^Xc\yc)7b^"'\b^^\bi^^^ and (^Xc\yc),b''''^\b^i3\bif^^^, where the indices, a 
and s, denote absolute value and square, respectively. 

The next iteration is in connection with a 5-dimension hyperrectangle, 
which has the following features: (1) it is centered near the previously de- 
termined point of minimum; (2) it is reduced in size by a factor of about 10, 
provided inequalities (79), (80), and (81) continue to hold; and (3) there is 
no point in common with the earlier hypergrid. Then two additional points 
of minimum are calculated and the next iteration is allowed to start. 

The computations end when the sum of absolute values and squares of 
absolute logarithmic residuals fall below an assigned treshold, which yields: 

^|i?,|<e(") ; (82) 
1=1 

j:Rl<e^'^ ; (83) 

i=l 

where the sum is performed on the range of interest, defined by Eq. (67). 
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